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ABSTRACT 

We present measurements of the Ca n triplet and the Ca4227 Lick-index for a sample of early- type 
galaxies in the Coma cluster, deriving, for the first time, their corresponding relationships with velocity 
dispersion. Compared with a similar subsample of ellipticals in the field, Coma galaxies with velocity 
dispersions in the range w 180 — 270 km s^ 1 exhibit significant differences in the strengths of the Ca 
features, suggesting an influence of the environment on the star formation histories of these galaxies. We 
argue that the main scenarios previously proposed to explain the relatively low Can triplet of galaxies 
are not able by themselves to simultaneously reconcile the strengths of the two Ca indices in both 
environments. 

Subject headings: galaxies: evolution — galaxies: clusters — galaxies: formation — galaxies: 
abundances — galaxies: stellar content 
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1. INTRODUCTION 

The near-infrared (nIR) Can triplet constitutes the 
strongest absorption feature in the integrated spectra 
around 8600 A of relatively old stellar populations. Given 
its well-known capability to derive the metallicity of galac- 
tic globular clusters (Armandroff & Zinn 1988, hereafter 
A88; Rutledge, Hesser & Stetson 1997), this spectral fea- 
ture was considered a useful tool to study the stellar popu- 
lations of normal galaxies. However, the first observational 
work (Cohen 1979; Bica & Alloin 1987; Terlevich, Diaz 
& Terlevich 1990; Houdashelt 1995) reported only small 
differences among the integrated Call triplet of several 
types of galaxies, remarkably at variance with the trends 
found when using other metallicity indicators (e.g. the Mg 
triplet at 5175A). The difficulties of previous single stellar 
population (SSP) models to make reliable predictions of 
the Ca II triplet (due to severe deficiencies in the available 
stellar libraries) prevented an understanding of the above 
apparent inconsistency. 

After a series of papers devoted to construct realis- 
tic models for disentangling the actual behaviour of the 
Can triplet in stellar populations (Cenarro et al. 2001a, 
hereafter C01; Cenarro et al. 2001b; Cenarro et al. 2002, 
hereafter C02; Vazdekis et al. 2003, hereafter V03), the 
above topic is being revisited and constitutes a controver- 
sial matter of debate. There exist recent evidences showing 
that, unlike other metal-lines, the Can triplet is surpris- 
ingly anti-correlated with central velocity dispersion (ctq) 
for several types of galaxies: ellipticals (Es) in the field 
(Saglia et al. 2002, hereafter S02; Cenarro et al. 2003, here- 
after C03), bulges of spirals (Falcon-Barroso et al. 2003, 
hereafter F03) and dwarf Es (dEs) in the Fornax cluster 
(Michielsen et al. 2003). Even more importantly, the mea- 
sured values in giant Es lie well below model predictions 
for any reasonable choice of stellar populations parameters 
(age and metallicity). A wide range of interpretations, like 
the existence of a dwarf-enriched stellar population, Ca 



underabundances, Ca depletion in the interstellar medium 
or a composite stellar population have been discussed in 
the above references. Unfortunately, the lack of definitive 
evidences makes no consensus to exist at the moment. 

In the blue spectral range, Ca4227 is considered as the 
most Ca-sensitive Lick-index (Tripicco & Bell 1995). In- 
triguingly, despite Ca is an a-element like Mg, the metal- 
licities inferred for Es using Ca4227 suggest that Ca is 
not enhanced -even it could be depressed- with respect 
to Fe (Worthey 1992; Vazdekis et al. 1997, hereafter V97; 
Worthey 1998; Trager et al. 1998; Vazdekis et al. 2001; 
Proctor & Sansom 2002, hereafter P02). Once again, the 
existence of Ca underabundances have been demanded to 
explain the low Ca4227 values (Thomas, Maraston & Ben- 
der 2003b; hereafter T03b). 

It is well known that the environment is expected to 
play a decisive role in the assembling and star formation 
history of galaxies. Hierarchical scenarios predict that 
Es in rich clusters were mostly formed at high redshifts, 
whereas field Es may have experienced an extended and 
more complex star formation history (Kauffmann & Char- 
lot 1998). Although no differences in the [Mg/Fe] ratios 
seem to exist between cluster and field Es (J0rgensen 1999; 
Kuntschner et al. 2002), Sanchez-Blazquez et al. (2003, 
hereafter S03) have recently found significant differences 
in C and N abundance ratios between Coma and field Es. 
These differences impose strong constraints to models of 
galaxy formation and chemical evolution. 

The above picture motivated us to extend the previous 
work to the analysis of the Can triplet and the Ca4227 
index in galaxies of different environments. First of all, 
because it might reveal new clues on the so confusing 
Can triplet behaviour. Secondly, because the fact that 
Ca seems to be an anomalous a-element makes it specially 
interesting from the point of view of chemical evolution. 

2. OBSERVATIONS AND DATA REDUCTION 
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This work makes use of long-slit spectroscopic data for 
two samples of galaxies: Es from the Coma cluster and 
Es in the field. Both were observed during two different 
runs (1999 March and 2001 April) using ISIS, the double- 
arm spectrograph at the 4.2 m William Herschcl Telescope 
(Observatorio del Roque de los Muchachos, La Palma). 
All details about the sample, observations and data re- 
duction for Es in the field are given in C03 and Sanchez- 
Blazqucz (2004; hereafter S04), for the nIR and blue-arm 
data respectively. We also refer the reader to S04 for a 
description of the reduction of the blue spectroscopic data 
for the Coma galaxies. In this section we just concentrate 
on their nIR spectral range (8355 - 9164 A). 

Our sample of galaxies from the Coma cluster consists 
of 28 early-type galaxies (E - SO) spanning a wide range of 
central velocity dispersions (80 < <r < 370kms _1 ). We 
have rejected the faintest low-mass galaxies in S03 because 
of quite large errors. When it was feasible, the 3.5arcmin 
long slit (2arcsec width, providing 2.9 A FWHM spec- 
tral resolution) was rotated to include two galaxies in the 
same exposure. Otherwise, single galaxies were observed 
by aligning the slit with the major axis. Exposure times 
of 1200 - 2000 s per galaxy leaded to signal-to-noise ratios 
per angstrom from 30 to 120 A" 1 in a central aperture of 
radius i? ff/8. 

Standard spectroscopic reduction procedures were per- 
formed with REDmE (Cardicl 1999). Apart from an ac- 
curate correction for the fringe pattern, we took special 
care on the sky subtraction and the correction for telluric 
absorptions at A > 8950 A (by using the telluric pattern 
derived from the spectra of flux standard stars). Both 
are specially critical for Coma galaxies since the redshift 
of the cluster (z ~ 0.025) places some index bandpasses in 
regions affected by atmospheric absorptions and strong sky 
emission lines. The availability of error spectra for each 
galaxy frame allowed us to estimate reliable uncertainties 
in the measurements of the indices. Index errors account 
for photon statistics, uncertainties in the flux calibration 
and radial velocity determinations (see details in C01). 

Final spectra were relative-flux calibrated using spec- 
trophotometric standard stars observed several times at 
different air masses. In order to transform the nIR spec- 
tra to the V03 spectrophotometric system, a sample of 39 
stars (from B to late M spectral types) in common with 
C01 were observed during twilights. Following C03, they 
were also employed as templates for <ro determinations. 
Rather than using template-dependent polynomials, sys- 
tematic differences between indices measured at different 
spectral resolutions were prevented by broadening all the 
spectra up to the largest Co of the galaxy sample, namely 
370 km s^ 1 , corresponding to NGC4889. 

The observation of four Es in common with C03 allowed 
us to confirm an homogeneous spectrophotometric system 
for the blue and nIR spectra of both observing runs. 

3. INDEX— (To RELATIONS 

We analyze, for the first time, the behaviours of the Ca II 
triplet, the H Paschen series (CaT*, CaT and PaT line- 
strength indices; sec definitions in C01) and the Ca4227 
Lick-index as a function of Co for Es in the Coma cluster. 
All the measurements correspond to a central aperture of 

3 http : //www . ucm . es/inf o/Astrof /ellipt /CATRIPLET . html 



radius R e g /8 (or 1 arcsec for galaxies with R c s < 8 arcsec) . 
There exists a database 3 listing the galaxy sample as well 
as the indices and cto determinations. 

Figure 1 shows the measurements of the above indices 
versus do for Es in the Coma cluster and, for comparisons, 
for a similar subsamplc of Es in the field. Error-weighted, 
least-squares linear fits to all data have been computed 
for both samples (see labels in Fig. 1). CaT* and CaT in 
Coma Es exhibit an anti-correlation with ctq (panels a and 
b) which is compatible with the one derived for field Es 
(C03). However, because of the errors of the computed 
slopes, we cannot rule out that a flat trend with a could 
exist for the Coma sample. As in the case of field Es, a flat 
behaviour is derived for the PaT index (panel c), whilst 
both slightly positive and negative (nearly flat) trends are 
evidenced for the Ca4227 indices of field and Coma Es 
respectively (panel d). In brief, as regards to the qualita- 
tive behaviours with o~o of the above indices, we find no 
important differences between field and Coma Es. 

The most striking differences between both samples 
arise from the mean absolute values of their Ca indices. 
In particular, though intermediate-low mass Es (er < 
180 km s -1 ) in both samples exhibit similar mean index 
values, we find clear differences for those Es within the 
range 180 < gq < 270 kms -1 (in the following we will 
refer to them as massive Es, the ones being the main sub- 
ject of discussion in this letter). Note that, other than 
visual inspection, there exist not any a priori argument to 
justify the above subsample selection. Interestingly, the 
fact that the differences in C4668 and CN 2 for Coma and 
field Es are also particularly evident in a similar range of 
velocity dispersion (S03) supports the idea that massive 
Es in different environments could indeed have different 
stellar population properties. Unfortunatelly, the scarcity 
of Es with cto > 270 kms -1 in our Coma sample prevents 
us to confirm whether the above behaviour extends to the 
high-mass end of the Es family. 

For this subsample of massive galaxies consisting of 17 
Coma Es and 15 Es in the field, we have computed error- 
weighted means of the indices and their corresponding 
standard errors (see labels in Fig. 1). The statistical signif- 
icance of the differences between the mean indices of both 
samples is studied by using a t-test (see Fig. 1; for a signif- 
icance level of a = 0.05, t values larger than 1.96 indicate 
that a significant difference exists). Apart from the fact 
that no significant differences are found for the PaT in- 
dex (as it is expected for relatively old stellar populations; 
see V03), we conclude the existence of a intriguing be- 
haviour: while field Es exhibit Ca4227 values significantly 
larger than Coma Es, the contrary occurs for the CaT* 
and CaT indices. In order to check that the above results 
are not driven by uncertainties in the correction for nIR 
telluric absorptions, we also measured the Ca2 and Ca3 
classical indices by A88 (Ca2+3 = Ca2 + Ca3). The high 
sensitivities of these indices to velocity dispersion and to 
the presence of Paschen lines (see C01) are not important 
when, as it is the case, they are measured on the spectra 
of old stellar populations at equal spectral resolution. In 
turn, the major improvement lies in the fact that their 
continuum sidebands are located in regions free from tel- 
luric absorptions for the redshift of the cluster. The data 
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in Fig. 1 (panel e) prove that the trends and differences 
inferred from these indices are compatible with the ones 
derived for CaT* and CaT. 

4. DISCUSSION 

As it has been recently reported by S02, C03 and F03, 
single and composite stellar populations models with so- 
lar abundance ratios ([a/Fe]= 0.0) and a standard initial 
mass function (IMF) cannot reproduce the low Call values 
of massive galaxies by solely assuming age and metallicity 
variations. Among others, two alternate scenarios have 
been proposed: i) the existence of Ca underabundances 
and ii) the existence of a dwarf-enriched stellar population. 
In the following paragraphs we explore in turn whether the 
above scenarios are consistent enough to explain simulta- 
neously the discrepancy with the model predictions and 
the differences between distinct environments. 

i) Even though there exist no evolutionary synthesis 
models for the Ca II triplet in which non-solar abundances 
ratios are considered, the Ca underabundant hypothesis 
has been proposed as a plausible scenario to explain the 
low CaT of massive galaxies (S02; F03). Making use of 
stellar population models with variable element abundance 
ratios for the Lick indices (Thomas, Maraston & Ben- 
der 2003a; hereafter T03a), T03b suggest the existence 
of a Ca depletion, with respect to the rest of a— elements, 
that increases with galaxy mass. Under this scenario, and 
assuming that the integrated Ca n triplet indeed traces the 
Ca abundance, the relative differences in the CaT of mas- 
sive Coma and field Es should be naturally explained as 
differences in the levels of Ca underabundances, being field 
Es "more Ca underabundant" than Coma galaxies. How- 
ever, there exist a strong contradiction in the fact that 
while mean Can indices of massive Coma Es are larger 
than of those in the field, the contrary occurs for Ca4227. 
Furthermore, taking into account that Ca4227 is affected 
by C and N abundances in the sense that the stronger CN 
the lower Ca4227 (Schiavon et al. 2002; P02; T03a), and 
given that field Es exhibit CN overabundances larger than 
those in Coma Es (S03), the relative differences in Ca4227 
cannot be explained by just CN effects either. In fact, if we 
consider such an effect, the actual difference between the 
Ca4227 strengths of both samples should be even larger 
than that in Fig. 1 (panel d). In summary, just the exis- 
tence of Ca underabundances and/or CN effects is unable 
to reconcile Ca4227 and CaT* at the same time. 

ii) Since the Can triplet of giant stars is stronger than 
that of dwarf stars, the existence of a dwarf-enriched stel- 
lar population could account for the low values of massive 
Es. This scenario can be parametrized by means of a vari- 
ation of the IMF-slope (/i; assuming a Salpeter-likc IMF). 
The evolutionary synthesis models by V03 predict that, 
for old SSPs with intcrmcdiate-to-high metallicities, CaT* 
and CaT mainly depend on /i (they decrease with the in- 
creasing fx because of a larger dwarf-to-giant stars ratio) . 
In the light of the above predictions, massive Es are inter- 
preted to have both metallicities and dwarf-to-giant ratios 
larger than low mass Es (C03). 

In Figure 2 we make use of the models by Vazdekis 
(1999; hereafter V99) and V03 to compare Coma and field 
Es under this scenario. The only way to reconcile H/3, 
[MgFe]' (the one being independent of typical a /Fe ratios 



as defined by T03a) and CaT* at the same time is invok- 
ing different ages, metallicities and dwarf-to-giant ratios. 
In turn, the later also affects the derived ages and metal- 
licities because of the non-negligible dependendence of H/3 
to the dwarf-to-giant ratio. From panels a, b, c and d, 
we derive similar metallicities (above solar) for both sub- 
samples, a mean age for field Es lower (~ 8.9 Gyr) than 
that for Coma Es (~ 11.2 Gyr), and mean dwarf-to-giant 
ratios larger for the former (/i ~ 2.8) than for the latter 
(fi ~ 1.8). However, when this scenario is simultaneously 
applied to Ca4227, the model predictions fail to reproduce 
the locus of the galaxies (panel c). Therefore, even though 
relative differences between the Ca indices of Coma and 
field Es could be compatible with just variations of the 
dwarf-to-giant ratio (the larger /x, the lower CaT* and the 
larger Ca4227; panels c, d, e and /), this argument is not 
enough to reconcile absolute values of CaT* and Ca4227 
by itself. 

At this point, one notices that the current status is 
not quite optimistic. Could it be that both scenarios co- 
exist driving an apparent inconsistency? Even though 
the sensitivity of Ca4227 to unusual dwarf-to-giant ra- 
tios is not negligible, age and metal abundances effects 
are mainly expected to govern its behaviour. The arrows 
in panels e and / of Fig. 2 would indicate the offsets in 
Ca4227 -arising not only from the existence of Ca sub- 
abundances but also from possible dilution effects due to 
CN enhancements- required to reproduce the SSPs prop- 
erties in the previous panels. On the contrary, several 
arguments support the idea that the Ca II triplet is rather 
unsensitive to Ca abundances. First, the CaT index does 
not correlate with stellar [Ca/Fe] abundance ratios (C02). 
Also, for SSPs older than ~ 3 Gyr with intermediate-to- 
high metallicities, the integrated Can triplet saturates (or 
even decreases for dwarf-enriched SSPs) with the increas- 
ing metallicity (V03). This is indeed a by-product of the 
increasing luminosity-weight of late-M stellar types (with 
very weak Ca n lines and strong TiO bands in their spec- 
tra) as the temperature of the isochrone becomes colder. 
Therefore, just because of temperature effects, very metal 
rich stellar populations with high Ca abundances are ex- 
pected to exhibit CaT values lower than if, for example, 
they were solar-abundant (see more details in V03). In 
this sense, Can triplet indices are expected to be sensitive 
to global metallicity rather than to Ca abundances. 

To conclude, we report, for the first time, systematic dif- 
ferences between the integrated Can triplet and Ca4227 
of Es in the Coma cluster and in the field which are dif- 
ficult to reconcile at the same time. Different abundance 
ratios and/or dwarf-enriched stellar populations could be 
conspiring to produce apparently inconsistent Ca indices. 
As proposed in S03, these kind of differences are proba- 
bly driven by the effect of the environment on chemical 
evolution and galaxy formation. Unfortunately, far from 
concluding any definitive prove for distinct star formation 
histories in massive Es located at different environments, 
the fact is that the current picture of calcium in galaxies 
seems to go beyond the capabilities of the available SSPs 
models. An extensive effort in the field of nucleosynthesis 
and stellar atmospheres are demanded to understand how 
different abundance ratios affect, not only line-strength in- 
dices, but also the physics and temperature of isochrones. 
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In the mean time, the collection of more and high-quality 
data of Ca line-strengths is crucial for trying to understand 
the controversial Ca puzzle. 
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Fig. 1. — CaT*, CaT, PaT, Ca4227 and Ca2+3(= Ca2 + Ca3 indices by A88) versus central velocity dispersion for Es in the Coma 
cluster (triangles) and in the field (circles). For comparisons, the indices of 4 galaxies in common between both runs are included (open 
circles and filled triangle). All the indices are measured at 370 kms -1 spectral resolution and corrected to the system defined by the models 
(V03). Solid (Coma) and dotted (field) lines represent error-weighted, least-squares linear fits to all data (7 = A + B ao; see the labels, 
with C and F subscripts refering to Coma and field Es). Iq and 7p are error- weighted mean indices for Coma and field Es within the range 
180 < cro < 270 kms -1 (shaded region). A t value larger than ~ 1.96 indicates that Iq and Ip are statistically different. NGC4742 (see arrow 
in panel d), with a very low Ca4227 because of a central young stellar population (see C03), has been rejected from the linear fits and the 
mean index computations. 
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Fig. 2. — Mean values of H/3, [MgFe]', Ca4227 and CaT* for Coma (open triangle) and field (filled circle) Es within the range 180 < <tq < 
270 kms -1 . Rather than typical errors, error bars correspond to the r.m.s. standard deviation of the data (larger than the former ones). SSPs 
model predictions by V03 and V99 at 370kms~ 1 spectral resolution respectively are overplotted. In panels a and b, age and metallicity vary 
from 5.62 to 14.13 Gyr and from —0.68 to +0.20 dex at different Salpctcr-likc IMF slopes (fi; see the labels). Panels c, d, e and / exhibit 
model predictions at fixed age (see the labels) with varying fi (0.3 — 3.3) and metallicity (as in panels a and b). Arrows in panels e and / are 
explained in the text. 



